Demonstration of Quadrature Squeezed Surface-Plasmons in a Gold Waveguide 
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We report on the efficient generation, propagation, and re-emission of squeezed long-range surface- 
plasmon polaritons (SPPs) in a gold waveguide. Squeezed light is used to excite the non-classical 
SPPs and the re-emitted quantum state is fully quantum characterized by complete tomographic 
reconstruction of the density matrix. We find that the plasmon-assisted transmission of non-classical 
light in metallic waveguides can be described by a Hamiltonian analogue to a beam splitter. This 
result is explained theoretically. 

PACS numbers: 03.67.-a, 42.50.Lc, 42.50.Nn, 73.20.Mf 



Enormous interest has recently been devoted to the 
emergent field of quantum plasmonics due to its unique 
capabilities in the way electromagnetic radiation can be 
localized and manipulated at the nanoscale. In partic- 
ular, integrated quantum technologies based on surface 
plasmons hold great promises for quantum information 
processing, since it allows for scalability, miniaturization, 
and coherent coupling to single emitters P, |3, S, 0, 0] ■ 
To enable these quantum information processing tech- 
nologies with high fidelity, it is of paramount impor- 
tance, that the nonclassicality of the plasmonic modes 
is preserved in propagation. The first experiment ver- 
ifying the preservation of entanglement in plasmonic 
nanostructures was carried out by Alterwischer et al. Q . 
They demonstrated the survival of polarization entan- 
glement after plasmonic propagation through subwave- 
length holes in a metal film. The preservation of energy 
time entanglement in a perforated metal film as well as 
in a thin conducting waveguide was later demonstrated 
by Fasel et al. ■ These experiments have witnessed the 
preservation of probabilistically prepared entanglement 
(thus neglecting the vacuum contribution) described in 
the two dimensional Hilbert space. 

In the present Letter we investigate the compatibil- 
ity of the quantum plasmonic technology with the con- 
tinuous variable quantum domain (described in the in- 
finite dimensional Hilbert space) by demonstrating the 
plasmonic excitation, propagation, and detection of de- 
terministically prepared quadrature squeezed vacuum 
states. We show that a squeezed vacuum state excite 
an electron resonance on the surfaces of a metallic gold 
waveguide to form a surface plasmon polariton (SPP). 
Despite loss and decoherence in the plasmonic mode we 
demonstrate that quadrature squeezing is retained in the 
retrieved light state. Importantly, we fully characterize 
the input state and output state by performing a com- 
plete quantum tomographic reconstruction of the states 



density matrix. This is in strong contrast to previous 
experiments on plasmon assisted quantum state trans- 
mission, where only a certain property of the quantum 
state was investigated. 

SPPs are combined electron oscillations and electro- 
magnetic waves propagating along the interface between 
a conductor and a dielectric medium 0]. By reducing 
the thickness of the metal film the SPP from the upper 
interface and the lower interface can couple. This cou- 
pling results in the formation of either long-range (LR) 
SPPs or short-range (SR) SPPs depending on the propa- 
gation length of the mode Q. In general, the coupled 
LR-SPP (SR-SPP) modes are characterized by an in- 
creased (decreased) propagation length compared to a 
SPP propagating along a single interface. LR-SPPs (SR- 
SPPs) are further characterized by a symmetric (anti- 
symmetric) longitudinal electric field distribution, which 
allows the LR-SPPs to be efficiently excited using end- 
fire coupling 0, . 

A sketch of the sample as well as the experimental 
setup is shown in Fig.[lja). The sample consists of a gold 
stripe embedded in lossless transparent polymer Benzo- 
cyclobutene (BCB). It was prepared by first spinning a 
12-15/iTO layer of BCB onto a silicon substrate. The 
metallic waveguide was then patterned using standard 
UV lithography, followed by gold deposition by means of 
electron beam evaporation and lift-off. Finally, another 
layer of BCB with a thickness of 10 /im was spun onto 
the sample with exactly the same spinning conditions 
to ensure a symmetric dielectric environment of the gold 
stripe. In Fig.[ljb) we present a power density plot of the 
LR-SPP mode in the transverse plane of the gold stripe. 
This was obtained by a finite element method simulation 
for the metal stripe geometry as presented in Fig. [ija). 
The arrows in Fig. \V[h) show the direction of the electric 
field vector in the transverse plane of the sample, which 
is, with respect to the metal stripe surface, transverse 



2 



magnetic (TM) polarized. This fact is in agreement with 
the analytical model for a sample of infinite width Q . 
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FIG. 1: (Color online), (a) Sketch of the experimental setup 
and LR-SPP sample: The sample is made of gold embedded 
in lossless transparent BCB and its dimensions are specified 
as written in the figure. Experimental setup: OPO - optical 
parametric oscillator, PPKTP - periodically poled KTP crys- 
tal, SV - squeezed vacuum, A/2 - half- wave plate, $ - piezo 
actuated mirror for phase variation, 50:50 - symmetric beam 
splitter, and HD - homodyne detection scheme, (b) Normal- 
ized power density plot of the LR-SPP mode for sample di- 
mensions as specified in (a). The arrows show the direction 
of the electric field in the transverse plane. At the experi- 
mental wavelength of 1064nm the refractive indices of gold 
and BCB are specified to be uau = 0.2381 + i7.7199 [ij] and 
nBCB = 1.539, which were used in the simulation, (c) Depen- 
dence of the total transmission T on the angle a of the linear 
polarized incident light field, a is measured with respect to 
TM polarization and set by the A/2 wave plate. 

As a source of non-classical light we use a bow-tie 
shaped optical parametric oscillator (OPO) operating be- 
low threshold, see Fig.[T]Ja) [12]. The non-linear medium 
inside the OPO is a type-Ill periodically poled KTP (PP- 
KTP) crystal of 10mm length, which is placed at the 
beam waist between two curved highly reflective mirrors 
with 25mm radius of curvature. We pump the PPKTP 
crystal with a frequency doubled NdiYAG laser field hav- 



ing a wavelength of 532nm, thus producing squeezed light 
centered around 1064nm. The round trip path length of 
the OPO cavity is 275mm, it has a free spectral range 
of l.lGHz, and a bandwidth of approximately 21MHz. 
The out-coupling mirror has a power transmissivity of 
T=10%. 

Before investigating the quantum properties of a SPP, 
we characterize its classical properties. For this purpose, 
we seed the OPO cavity with an auxiliary light field at 
1064nm with the OPO pump beam blocked. A half wave 
plate is placed between the OPO and the gold sample to 
control the linear polarization of the light field. The light 
field is then focused onto the sample with an asphcric lens 
ensuring a good mode matching to the LR-SPP mode. 
At maximum, we have observed a total transmission of 
32.8 ±0.5% for TM polarized incident Hght. When rotat- 
ing the polarization of the incident light field, the trans- 
mission decreases as cos^(a), where a is the angle of the 
half-wave plate with respect to the TM direction. For 
a = 90° only scattered as well as higher order modes 
are measured. This result clearly demonstrates the ex- 
pected polarization dependence and thus proves, that the 
LR-SPP mode is excited 0. 

Having established the classical excitation of a SPP, we 
now proceed with the analysis of non-classical SPPs. In 
the following we fully characterize the quantum state of 
light before it is launched into the gold stripe, and also 
after transmission and re-emission. The results of the 
analysis are shown in Fig. [2l For the characterization of 
the input state, the sample is removed from the setup 
and the squeezed vacuum light field is measured with a 
homodyne detection system, as shown in Fig. [Ija). The 
difference signal of the two detectors is amplified, mixed 
down at an optical sideband frequency of 4.7MIIz, low 
pass filtered with a bandwidth of 150kHz, and finally dig- 
itized with 8bit resolution. A quadrature measurement 
set containing more than 5.5M data points is shown in 
Fig. [2l^a) for a rotation in phase space from to 27r. By 
applying the maximum likelihood algorithm 1J| on this 
data set we reconstruct the squeezed vacuum states den- 
sity matrix pi„, whose absolute values are presented in 
Fig. Wljo) ■ From the density matrix pi„ , the Wigner func- 
tion W{X, P) has been reconstructed, which is presented 
in Fig. [ll^c). For various qudratures, the noise power can 
either be calculated directly from the time resolved data, 
as shown by the black solid line in Fig.[2I^d), or from the 
states density matrix as shown by the red dashed line 
in Fig. [2Ud) . The degree of squeezing and anti-squeezing 
are found to be -1.9±0.1dB and 6.1±0.1dB, respectively. 

The squeezed light is then carefully injected into the 
gold stripe, thereby linearly mapping the non-classical 
photon statistics onto a plasmonic state, as described for 
the exciation with a single photon in Ref. [l^]- Subse- 
quently, the squeezed SPP propagates through the sam- 
ple, is efficiently re-emitted at the rear end, and mea- 
sured with the homodyne detector. The time resolved 
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FIG. 2: (Color online). Experimental results of the squeezed 
vacuum state (input) and the quantum state after LR-SPP 
excitation, propagation, and reemission (output), (a) and (e) 
Time resolved data for various electric field quadratures while 
scanning the relative phase between LO and signal from to 
2n. The visibility between signal light field and LO is 88.5% 
in case of the squeezed vacuum field and 90.3% in case of the 
LR-SPP mode. The data are normalized in a basis, where 
[x,p] — i. (b) and (f) Reconstructed density matrices pi„ and 
Plr-spp in the photon number basis obtained by applying 
the maximum likelihood algorithm to the data presented in 
(a) and (e), respectively, p has been reconstructed for a trun- 
cated Hilbert space with a photon number of n=15. (c) and 
(g) Calculated Wigner functions W{X, P) for the density ma- 
trices shown in (b) and (f), respectively, (d) and (h) Noise 
power of the data set presented in (a) and (f ) (black solid line) 
and calculated from the density matrices (red dashed line). 
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FIG. 3: Fidelity between beam splitter mode poutijf) obtained 
after applying the beam splitter operator UBsiv) on the input 
mode pin and LR-SPP output mode Plr-spp versus trans- 
missions rj. 



quadrature data are recorded similarly as was described 
above and the result is shown in Fig. ^e). From these 
data we reconstruct the density matrix plr-spp a-nd the 
Wigner function W{X, P), as illustrated in Fig. [^ff) and 
(g), respectively. By calculating the quadrature variances 
(either directly from the time-rcsolvcd data or from the 
density matrix) we find a minimum of -0.7 ± O.ldB and 
a maximum of 3.2 ± O.ldB with respect to the shot noise 
level, as is presented in Fig.jSl^h). We therefore conclude, 
that the retrieved light state is squeezed and that the 
squeezing survived the plasmonic propagation. 

We furthermore investigate whether the operation, 
that transforms the density matrix of the input state 
Pin to that of the output state plr-spp can be de- 
scribed by the unitary beam splitter operator Ubs ~ 
exp{^{a) be^'^ — afe^e"**)}. Here, d and b are the field 
operators of the beam splitter input modes, $ is the rel- 
ative phase between the modes a and 6, and 6 is linked to 
the transmission rj via -q = cos(|). The expected output 

state is thus Poutiji) = Tr{UBs{v)Pin \Q){Q\Ulsiv)}: 
where the trace is taken over one of the output modes of 
the beam splitter. To see whether Poutii]) is similar to 
the actually measured output state Plr~spp we compute 
the fidelity J^(f?) between the two states. The fidelity is 
givenby f (77) = tr{ ^ pout {'n)PLR.-SPP\/ Pout (??)}^/^ [H, 
with < F(?7) < 1 and F{ri) = 1 if and only if 
Poutiv) = Plr-spp- In Fig. [3l wc plot the fidelity F{ri) 
versus the transmission rj of the beam splitter, which 
reaches F = 0.993 at maximum for 77 = 0.33. The large 
overlap between poutir] = 0.33) and plr-spp are in ex- 
cellent agreement with the measurements presented in 
Fig- dc) and proves, that the plasmonic decoherence is 
linear and thus can be simulated by a beam splitter in- 
teraction. 

In the following we justify this conclusion by using sim- 
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pie theoretical arguments. From the quantum mechan- 
ical Maxwell equations the propagation of the field in 
an (electrically) polarizable medium may in general be 
described by 



92 1 \ . 1 92 ^ 
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(1) 



where E is the electric field and P is the polariza- 
tion of the medium. We are here interested in the 
component at optical frequencies and split the electric 
field into the positive and negative frequency compo- 
nents E — -I- where the positive (negative) 
frequency component contains the photon annihilation 
(creation) operator. Similarly, we split the polarization 
P = p(+) -I- The positive frequency component of 
the polarization may be expanded in terms of the electric 
field 



(2) 



where we include the possibility that the expansion coeffi- 
cient Cran may depend on a set of state operators {Oi} de- 
scribing the state of the polarizable medium. During the 
propagation, the photonic excitations of the electric field 
are mapped onto the (quasiparticle) excitations of the po- 
larizable medium. The expansion coefficient Cmn{\Oi\) 
thus gives access to information about the nature of the 
excitation of the medium and, e.g., quasi particle interac- 
tions. Our experiment is performed in the low intensity 
limit, in which case the polarization reduces to p(+) = 
coo({'5i}) -I- Coi({Oi})i?^+\ where we have excluded the 
E^~^ term since a negligible intensity is observed when 
the incoming field is in vacuum. This also means that we 
must require (cJqCoo) = 0. Since the equation of motion 
^ is linear, the expression for P(+) means that the out- 
going field -Eq^j can be written as a combination of two 
terms ^^+' = G({Oz})coo({0,:}) + G({0.})4+\ where 
G({Oi}) is the Greens function of the plasmonic prop- 
agation. If we ignore any dependence on internal state 
operators {Oi} in the expansion coefficient cqo and cqi 
in Eq. ([2]), we find that the input/output relation for a 
single mode operator a is given by 



(3) 



where i) is a combination of the cqq operators with 
the factor -^/l — r\ separated out for convenience. From 
Coo) = then immediately follows that (wt?)) — 
and consistency of the commutation relations requires 



■-oo 



that 



= 1. w is thus a single mode vacuum op- 
erator in agreement with our experimental observation. 
Our results therefore show that the classical description 
p(+^ = (e — eo)^''"^'' (where the permittivity of the mate- 
rial e is just a constant) remains valid for the metal down 



to the level of single photons if we add a vacuum contribu- 
tion to the operator equations. This conclusion is consis- 
tent with the theoretical interpretation in Ref. of the 
experiments in Ref. 0], where it is concluded that the 
polarization degrees of freedom leaves no "which- way" 
information in the solid. The present work, however, ex- 
tends that conclusion by showing that also the presence 
or absence of a photon leaves no "which- way" informa- 
tion. 

In conclusion, we have experimentally demonstrated 
the efficient plasmon-assisted propagation of squeezed 
vacuum states in a gold waveguide. Through complete 
reconstruction of the density matrices of the input and 
output fields we found that the plasmonic propagation 
can be described by a unitary beam splitting operation. 
The squeezing was therefore coherently transferred from 
the light state to the plasmonic state and back to the light 
state, solely degraded by vacuum noise. This demon- 
strated robustness of continuous variable quantum states 
in plasmon propagation suggests that continuous vari- 
able quantum information processing based on surface 
plasmons is feasible. 
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